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Experience has shown that efficiency usually increases when sep- 
arate traffic systems are combined into a single system. For example, 
if Group A contains 10 trunks and Group B 8 trunks, there should be 
fewer blocked calls if A and B are combined into a single group of 18 
trunks. It is intuitively clear that the separate systems are less efficient 
because a call can be blocked in one when trunks are idle in the 
other. Teletraffic engineers and queuing theorists widely accept such 
efficiency principles and often assume that their mathematical proofs 
are either trivial or already in the literature. This is not the case for 
two fundamental problems that concern combining blocking systems 
(as in the example above) and combining delay systems. For the 
simplest models, each problem reduces to the proof of an inequality 
involving the corresponding classical Erlang function. Here the two 
inequalities are proved in two different ways by exploiting general 
stochastic comparison concepts: first, by monotone likelihood-ratio 
methods and, second, by sample-path or "coupling" methods. These 
methods not only yield the desired inequalities and stronger compar- 
isons for the simplest models, but also apply to general arrival 
processes and general service- time distributions. However, it is as- 
sumed that the service- time distributions are the same in the systems 
being combined. This common-distribution condition is crucial since 
it may be disadvantageous to combine systems with different service- 
time distributions. For instance, the adverse effect of infrequent long 
calls in one system on frequent short calls in the other system can 
outweigh the benefits of making the two groups of servers mutually 
accessible. 

I. INTRODUCTION AND SUMMARY 

From extensive experience in teletraffic engineering, it is well known 
that congestion can often be reduced by sharing resources. The block- 
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ing probability in a loss system and the average waiting time in a delay 
system are usually much less when separate facilities serving separate 
streams of traffic are combined to serve all the streams together. 
Alternatively, for a given level of congestion, fewer facilities are usually 
required to serve the streams together. Sometimes such results are 
trivial: Whenever the combined system may be managed as if it were 
in fact separate systems, the optimal performance of the combined 
system is at least as good as that of the separate systems. However, 
such management is not allowed in the models treated here. In any 
case, the efficiency of shared resources is certainly a fundamental 
principle of teletraffic engineering. 

The purpose of this paper is to establish versions of this efficiency 
principle mathematically. Our first two results verify conjectures by 
Arthurs and Stuck. 1 To state our first result, let L(s, A u) denote the 
stationary loss or overflow rate in an M/M/s loss system (no waiting 
room) with s servers, arrival rate A, and individual service rate a. (See 
Kleinrock 2 for background on the queuing models.) It is well known 
that L(s, A, u) = \B(s, a), where a = A/u and B(s, a) is the familiar 
Erlang blocking formula: 

B(s,a) = (a7s!)/i; <a*/k!); d) 

see Jagerman 1 and references there. The first efficiency principle we 
establish says that L (s, A, u) is a subadditive function of (s, A) for each 
fixed ju: 

Theorem 1: For all positive integers s\ and s 2 and all positive real 
numbers Ai, A2 and ju, 

L(si + s it A] + Ki, u) < L{s u Ai, u) + L{s 2 , A 2 , u). (2) 

This yields immediately that 



g(£s„ £<fc|s£ G( Si , ai ) 
V-i <=i / 1-] 



for each integer n, where G(s, a) = aB{s, a), which is the version of 
Theorem 1 actually conjectured by Arthurs and Stuck. 1 

Of course, Theorem 1 should not surprise teletraffic engineers, since 
it can be inferred from common tables and graphs, but it has apparently 
not been proved before. It appears that all previous mathematical 
results can be described as "one-parameter" results. The relation (2) 
has been deduced for special cases in which quantities such as the 
blocking probability or the load per server are held constant. For 
example, it is known that if one combines separate groups of j and k 
trunks, each operating at a blocking probability 0.01, then the new 
blocking probability will be less than 0.01 (or, alternatively, the com- 
bined system can handle an increased total load and retain the same 
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0.01 blocking probability); see p. 68 of Cooper." Such results are often 
presented without rigorous mathematical support. 

From Paul Burke we learned about another special case that has 
been known for a long time. It is not difficult to show that B (ts, ta) is 
strictly decreasing in t (see the appendix), from which (2) easily follows 
in the case Ai/S] = X 2 /s 2 . Herbert Shulman has also shown that 
Theorem 1 follows easily from the monotonicity of B (ts, ta) in t and 
the convexity of B (s, a) in s for s > 1, but such convexity has not yet 
been established (see the appendix). For further discussion of other 
related work, see Section 5.1 of Kleinrock/' 

To state our second result, let D(s, A, (i) denote the mean steady- 
state delay in an M/M/s queue with infinite waiting room, fcfs (first- 
come, first-served) queue discipline, s servers, arrival rate A, and 
individual service rate fi. It is well known that D(s, A, /i) = C(s, A//i)/ 
(sfi — A), where C(s, a) is the Erlang delay function: 

Cfc^- „ "V(s-D!(s -a) (3) 

J <a*/A!) + ay(* - 1)! <« - a) 

The following result establishes subadditivity of D(s, A, fi) as a function 
of (s, A) for each fixed fi. Note that 

[A,/(Ai + X 2 )]D(s u X h fi) + [A 2 /(A, + X 2 )]D(s 2 , A 2 , fi) 

is the overall average delay experienced in the separate systems 
because Ai/(Ai + A 2 ) is the long-run proportion of customers to enter 
the first system. 

Theorem 2: For all positive integers S\ and s 2 and all positive real 
numbers X h Aa, and fi, 

D(si + s 2 , Ai + A 2 , fi) 

< [A,/(A, + X 2 )]D(8 h A,, ft) + [A 2 /(A, + X 2 )]D(s 2 , A 2 , fi). (4) 

This yields immediately that 

\ n 

t) 



Ml s it i Oijal H(s„a,) 



for each integer n, where H(s, a) = aC(s, a)/(s — a), which is the 
version of Theorem 2 actually conjectured by Arthurs and Stuck. 1 

In order to prove Theorems 1 and 2, we found it convenient to prove 
stronger results. It is helpful to see how the loss rate L(s, A, /x) and the 
mean delay D(s, X, fi) are related to the steady-state number of 
customers in the system, say Q. In the M/M/s loss system 

L(s, A,ju) =X-(iEQ (5) 

because A is the arrival rate and fiEQ is the service-completion rate; 
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the loss rate is that fraction of the arrivals that are not served. In the 
M/M/s delay system 

\(D(s, A, u) + a" 1 ) = EQ (6) 

by virtue of the fundamental relation L = XW; see Stidham. 

Let Qi be the steady-state number of customers in the ith system 
(i = 1, 2) and let Q be the steady-state number of customers in the 
combined system. Then Theorem 1 is equivalent to 

EQ > EQx + EQ 2 (?) 

for the loss systems, and Theorem 2 is equivalent to 

EQ < EQx + EQ 2 (8) 

for the delay systems. 

Instead of comparing the means in (7) and (8), we prove Theorems 
1 and 2 by making more general stochastic comparisons. We do this in 
two different ways. Our first method of proof is to compare the 
distribution of Q with the distribution of Qi + Q 2 . It turns out to be 
very easy to establish an appropriate ordering for the entire distribu- 
tions, which in turn implies the desired inequality for the means. The 
appropriate order is the monotone likelihood-ratio ordering. We define 
this ordering and prove the more general theorems implying Theorems 
1 and 2 in Section II. 

Our second method of proof is to compare entire stochastic processes 
rather than just stationary distributions. As corollaries we obtain 
stochastic-order relations for the stationary distributions which in turn 
also imply the desired inequalities (7) and (8) for the means. This 
approach has the advantage that the arrival processes can be arbitrary 
rather than Poisson and the service-time distributions can be general 
instead of exponential. The argument is also remarkably simple. The 
idea in this approach is to construct artificially the two stochastic 
processes being compared on the same probability space. The construc- 
tion is carried out so that each stochastic process individually has the 
correct distribution (family of finite-dimensional distributions) as orig- 
inally specified. We choose a special joint distribution so that each 
sample path of one process always lies below the corresponding sample 
path of the other process. Because the construction is artificial, the 
joint distribution of the two processes is not directly meaningful, but 
it implies a strong stochastic ordering for the processes. Such special 
constructions have been used previously to compare queuing processes; 
see Sonderman, 7 Whitt, 8 Wolff, 9 and references there. In fact, the 
generalization of Theorem 2 is a direct consequence of Wolff's theorem 
and the other proofs involve similar reasoning. We present our results 
using this approach in Section III. 
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In Theorems 1 and 2 we assume equal service rates in the two 
systems. It is natural to ask whether extensions of (2) and (4) hold 
when the service rates are unequal. In Section IV we show that, with 
unequal service rates, combining resources need not be more efficient; 
in fact it can substantially degrade performance. Infrequent "bad" 
customers from one system can adversely affect a large number of 
"good" customers from the other system. 



II. MONOTONE LIKELIHOOD-RATIO COMPARISONS 

Let X and Y be random variables assuming values in the nonnegative 
integers. We say X is less than or equal to Y in the monotone 
likelihood-ratio ordering and write X < r Y if 

P(X«A + 1) P(Y = k + l) 



■< 



P{X = k) P(Y = k) 



(9) 



for all integers k\ see page 208 of Ferguson. 1 " We say X is stochastically 
less than or equal to Y and write X <„, Y if Ef(X) < Ef(Y) for all 
nondecreasing real-valued functions f for which the expectations are 
well defined. Obviously, EX < EY whenever X <„, Y. What is important 
for Theorems 1 and 2 is that X < r Y implies X <„, Y. This is well 
known and not difficult to show. In fact, the monotone likelihood-ratio 
ordering is equivalent to stochastic order for all conditional distribu- 
tions obtained by conditioning on subsets, i.e., E(f(X)\X e A) ^ 
E(f(Y) \ Y G A) for all subsets A and all nondecreasing real- valued 
functions /; this property is discussed in Whitt 1112 ; see Keilson and 
Sumita 13 for additional material. 

Returning to the notation of (7) and (8), we obtain the following 
results which imply Theorems 1 and 2. 
Theorem 3: For the M/M/s loss systems, 

Qi + Qz <r Q- 
Theorem 4: For the M/M/s delay systems, 

Q<rQl + Q 2 . 

Theorems 3 and 4 can each be proved by simple calculations since 
the stationary distributions are known and easy to work with. To 
illustrate, we do one proof. 

Direct Proof of Theorem 3: Let a, = A,//i, for i = 1, 2. Then, using 
convolution, we obtain for some constant C 

P(Qi + Q 2 = k + l) = C X aM/iM 

i',+i.,-A+l 
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' X a'r'aWdi - «W 



A + 1 1SI]^6 

J!+I 2 -*+l 



+ .I!El I aWr l mh-D\ 



ft T 1 QS^SS] 



i,+i~=ft+l 



\k+l } 



0£I 2 SS 2 



P(Q = A) ^ D 

It is also significant that both Theorems 3 and 4 can be viewed as 
trivial corollaries of a more general theorem. This more general theo- 
rem is especially useful for comparisons when the limiting distributions 
are not known. To state our general result, consider two stochastic 
processes on the integers, Yi{t) and Y 2 (t), that move only by jumps up 
or down in unit steps to one of the neighboring states. Let all the 
transitions be governed by birth-and-death rates, but in contrast to 
those in birth-and-death processes, these rates may depend on infor- 
mation other than the current state such as the history of the process 
or other relevant variables. Let A,(A, I t ) and puik, It) be the birth-and- 
death rates, respectively, for the ith process (i = 1, 2) in state k with 
additional information J, at time L By having transitions governed by 
birth-and-death rates, we mean that 

P{Yi(t + h) - k + 1 1 Yi(t) = K h) - hki(t, I,) + o(h), 
P(Yi(t + h) = k - 1 | YM) = Kit) = hfiM, I,) + o(h) t 
and 

P(Yi(t + h) ~ k | Y t (t) - k, It) - 1 - h[\,(t, h) + M*. W + o(h), 

where o(h) means a quantity that converges to zero after division by 
h as a _». o. Let Xt and X 2 be random variables with the limiting 
distributions of these two stochastic processes, which we assume exist 
as proper distributions. Here is our general monotone likelihood-ratio 
comparison result. 

Theorem 5: Consider the processes Yi(t) and Y 2 (t) defined above. 
Suppose there exist sequences of constants {a,-(ft))} and {/?,(&))} such 
that 

A, (A, I,) < adk), Mk, I,) > a 2 {k) f 
jtii(A,I,)>j8,<ft), ftaft £)£&<*), 
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for all A and / ( . If ai(k)/fc(k + 1) < a 2 (A)/j8 2 (A + 1) for all A, then 

X.\ ^H r X 2 . 

Corollary: If 



A,(A,J,)<A a (A,/I) 



and 



fLi(k,I t ) >/i 2 (A, Id 
for all k, I t , and I' t , then Xi < r X 2 . 

Proof of Theorem 5: Look at the stationary flow between states k and 
k + 1. The flow from A to A + 1 is less than or equal to P(Xi = k)ai(k) 
for process 1 and greater than or equal to P{X2 = A)a 2 (A) for process 
2. Similarly, the stationary flow from A + 1 to k is greater than or 
equal to P{X\ = k + l)/?i(A) in process 1 and less than or equal to 
P(X 2 = k + 1)/?2(A) in process 2. Since the stationary flow from k to 
k + 1 must equal the stationary flow in the opposite direction, 

P{Xi = k)ai(k) > P(X, = k + 1)0,<A + 1) 
and 

P(X 2 = k)a 2 {k) < P(X 2 = k + l)Mk + 1). 

Consequently, 

P(X, - Jfe + 1) _ q,{A) ^ a 2 (A) ^ P(X 2 = A + 1) 
P(Xi = A) ~ ft (A + 1) ~ 2 (A + 1) ~ P(.Y 2 = A) " 

We can now apply the corollary to Theorem 5 to prove Theorems 3 
and A, 

Second Proof of Theorem 3: Note that the processes depicting the 
number of customers being served satisfy the hypotheses of Theorem 

5. In the case of two separate facilities, the sum is not a birth-and- 
death process because the rates depend not only on the total number 
but how many are in the individual facilities. When k customers are 
present, the death (service) rates are identical, but the birth (arrival) 
rates can be higher in the combined system because if one of the 
separate facilities is full, then it cannot accept any more arrivals. 
Hence, the hypotheses of the corollary to Theorem 5 are 
satisfied. D 

Proof of Theorem 4: Again we apply the corollary to Theorem 5. The 
reasoning is similar except here when A customers are present, the 
birth (arrival) rates are always identical, but the death (service) rates 
can be less with the separate facilities because there can be idle servers 
in one facility while there are customers waiting in the other 
facility. □ 
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III. SAMPLE PATH COMPARISONS 

Let {X(t), t > 0} and { Y(t) t t > 0} be real-valued stochastic processes. 
We call a real-valued function f defined on the space of all sample 
paths of X[t) and Y(t) nondecreasing ii f({x(t), t == 0}) < fi{y(t), t > 
0}) for all sample paths {x(t), t >: 0} and { y(t), t > 0} such that x(t) 
< y{t) for all t > 0. We say the stochastic process {X(t), t > 0} is 
stochastically less than or equal to the stochastic process { Y(t), t > 0} 
and write {X{t), t > 0} <* { Y(t), t > 0} if f({X(t), t > 0}) <„ f({Y(t), 
£> 0}) for all nondecreasing real-valued functions /"defined on the 
sample paths of X(t) and Y{t). Clearly, stochastic order of the processes 
implies X(t) < 8( Y(t) for each t [just use the projection: f({x(u), u > 
0)) = x(t)], but it is much stronger, applying to many other nonde- 
creasing functionals. In fact, since the queuing processes have sample 
paths with left and right limits everywhere, stochastic order of the 
processes is equivalent to stochastic order for all finite-dimensional 
(joint) distributions; see Section 4 of Kamae, Krengel, and O'Brien. 14 
Moreover, stochastic order of the processes here is equivalent to the 
possibility of a strong sample-path comparison. In particular, 

(X(*), * 2*0} =S* {Y(#),*&0} 

holds if and only if it is possible to construct stochastic processes 
{X(t), t > 0} and {Y(t), t > 0} on a common probability space such 
that {X(t), t > 0} has the same distribution as {X{t), t > 0}, 
{Y(t), t S 0} has the same distribution as {Y(t), t > 0}, and every 
sample path of {X(t) $ t > 0} lies below the corresponding sample path 
of {Y{t), t > 0}; see Theorem 1 of Kamae, Krengel, and O'Brien. 14 
What we do is apply the easy half of this equivalence— the fact that 
the sample-path construction implies stochastic order— to make sto- 
chastic comparisons between the queuing processes. The proofs here 
are done by actually constructing processes with the sample-path 
ordering. Previous uses of such constructions appear in Sonderman, 7 
Whitt, 8 Wolff, 9 and references therein. The approach is also closely 
related to the so-called "coupling" techniques; see Lindvall 15 and 
references therein. 

We begin with the generalization of Theorem 2 for delay systems 
because it follows directly from Wolff. 9 As before, we assume the fcfs 
discipline, but now we allow the arrival streams in the two separate 
systems to be arbitrary. We assume the service times are independent 
of the arrival processes and mutually independent and identically 
distributed, but they need not be exponentially distributed. Since the 
arrival process is assumed to be independent of the service-time 
sequence, the evolution of the arrival process cannot depend on the 
state of the system. This excludes finite-source models, for which 
counterexamples to the efficiency of sharing are easy to construct; for 
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example, see page 1377 of Benes. 16 Let Qi(t) be the number of customers 
in the ith system and let Q(t) be the number of customers in the 
combined system at time t. 
Theorem 6: (Wolff) If Qi(0) = Q 2 (0) = Q(0) = 0, then {Q(t), t > 0} 

Remarks: (i) Wolff 9 was actually interested in comparing the fcfs 
discipline with the cyclic assignment discipline in a single delay-system. 
He showed that the queue length process with the fcfs discipline is 
stochastically less than the queue length process in the same system 
with any other discipline. This result applies here because the two 
separate facilities can be interpreted as a single system with a special 
queue discipline: Just label the arrivals in the special system according 
to the stream from which they came and then assign them according 
to the fcfs discipline to one of the servers in the corresponding 
subgroup of servers. 

(ii) We can obtain corresponding results if the systems are not 
empty initially. For more general initial conditions, we can assume 
appropriate stochastic order for the residual service times at t = 0. 

(Hi) Wolff 9 also obtained similar comparison results for other pro- 
cesses, all of which hold here too: the departure epochs, the number of 
customers in queue, the total work (in service time) in the system, and 
the total work in queue. By the sample-path construction, the sto- 
chastic order jointly holds for all these processes. See Theorem 8 here. 

(iv) As a consequence of Theorem 6, Q(t) < s/ Q\{t) + Q2U) for each 
t. With the general assumptions here, steady-state distributions need 
not exist, but if Qi(t) and Q(t) converge in law to Q, and Q, respectively, 
as t — ► oo, then Q < s( Qi + Q 2 ; see Proposition 3 of Kamae, Krengel, 
and O'Brien. 14 The convergence of course holds in the setting of 
Theorem 2, so Theorem 6 implies (8) and thus Theorem 2. 

(v) Since Theorem 2 concerns the mean- waiting time, it is natural 
to ask if the steady-state waiting-time distribution is also stochastically 
less in the combined system. Unfortunately, in general it is not. The 
counterexample in Whitt 17 applies here too; the cyclic discipline there 
can be interpreted as arrivals to separate facilities. 

(vi) When the arrival streams are not Poisson, which we now permit, 
a new phenomenon occurs. Then the customers in the different streams 
experience different congestion when the systems are combined, even 
if the service times are independent and identically distributed. This 
phenomenon can be an important consideration in combining systems, 
but we do not consider it here; it has been studied by Kuczura. lH,m 

We now turn to our generalization of Theorem 1 for loss systems. In 
addition to allowing arbitrary arrival streams and general service-time 
distributions, we allow a finite waiting room. The number of waiting 
spaces in the combined system is the sum of the numbers of waiting 
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spaces in the separate systems. Let Ni(t) [N(t)] be the number of 
customers lost in the interval (0, t) in the ith separate system (in the 
combined system); let SM) [S(t)] be the number of service completions 
in the interval (0, t) in the ith separate system {in the combined 
system); and let Ct(t) [C(t)] be the amount of work performed— service 
given— in the interval (0, t) in the ith separate system {in the combined 
system). 

Theorem 7: If (?i(0) = Q 2 (0) = Q(0) = in these systems with finite 
waiting rooms, then 

{N(t), t > 0} <„, {NAt) + N 2 (t), t => 0}, 

{S(t) t t > 0} > K , {S,U) + S-At), t > 0}, 

and 

{C(t), t > 0} > s , {CM + Cs(t), t > 0}. 

Now assume that Ni{t)/t and N(t)/t converge {either in probability 
or with probability one) as t -+ oo. Let the limits be denoted L(s it k it 
Ai(t), F) and L(si + s 2 , *i + k 2 , A^t) + A 2 (t), F), respectively, with ki 
denoting the number of waiting spaces, Ai(t) the arbitrary arrival 
process and F(x) the general service-time c.d.f. From Theorem 7 we 
immediately obtain the following generalization of Theorem 1. 
Corollary: For all positive integers s u s 2 , ki and A 2 ; all arrival 
processes Ai(t) and A 2 (t); and all service time c.d.f.'s F(x) such that 
the loss-rate limits exist, 

L{si + s 2 , ki + k 2 , Ai(t) + A 2 (t), F) 

< L(si, k u AAt), F) + L(s 2 , k 2 , A 2 (t), F). 

To prove Theorem 7, we establish a finite-waiting-room generaliza- 
tion of Wolff's 9 comparison theorem. Following Wolff, we shall state 
the result in terms of the sample-path comparison. Since the joint 
distribution of the two systems being compared is artificially obtained, 
the appropriate conclusion is the general stochastic order as in Theo- 
rems 6 and 7. 

We carry out the artificial construction by letting the systems being 
compared have identical arrival processes and service times. Note that 
we are now focusing on a single (arbitrary) sample path. We let the 
nth service time v n be associated with the nth customer to enter 
service in each system rather than the nth arrival. Let a n be the arrival 
epoch of the nth arrival, < d < a 2 < • • • . We assume there are s 
servers operating in parallel and k extra waiting spaces in both systems. 
We also assume the systems are initially empty. 

One system, called the original system, will be the conventional 
system where the servers are fed by a single queue using a fcfs 
discipline. Moreover, there are k extra waiting spaces and arriving 
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customers enter the system if the number of customers in the system 
is less than s + k, and are lost otherwise. 

The other system, called the modified system, is any alternative to 
the original system which assigns customers to servers in some manner, 
independent of the sequence of service times {v„}, and which loses 
arrivals whenever the system is full and in some manner otherwise. 

Let a n be the arrival epoch of the rath customer. For the original 
system, let t„ be the time that the nth customer to enter the system 
enters; obviously t n = a )t for some k, k > n. Also, for the original 
system, let b n be the time the nth customer to begin service begins and 
let d n be the nth ordered departure epoch from the system. Let a n , t' n , 
bn, and d n be the corresponding quantities for the modified system. 
Theorem 8: For all integers n, t n < t n , b n < b' n , and d n < d' n . 
Proof: The sets of unordered departure epochs in the two systems are 
clearly {{b n + v n )} and {(b' n + v n )}, respectively. For the original 
system, 

dx = min {b, + u,} = min {£, + u,}, 

d„ = nth-order statistic from {(b, + v<): i < n + s} (10) 

and 

b n = max{*„, £*„_*}, n > 1, (11) 

where dj = if / = 0. For the modified system, 

dn = nth-order statistic from {(6J + vt): i < n + s) (12) 

and 

K > max{*„, dn-,), n > 1, (13) 

because in the modified system it is possible to have a positive queue 
and an idle server. 

Since the nth-order statistic is a monotonic function, to prove 
Theorem 8 it suffices to show that t n < t n and b n < 6J, for all n. We 
show this induction. Obviously bi = U = a, < t J < b[, 1 < ( < s. Suppose 
U < t\ and bi < b\ for all i, i £ n - 1. We first show that fe, < t' n . 
Suppose not; then 

tn>t' n > t'n-l > t n - h 

and thus n — 1 customers have entered both systems before the arrival 
associated with t' n . (Note that customers could arrive in batches, i.e., 
a k = a k +\ is a possibility, but this presents no serious difficulty.) 
However, by the induction hypothesis bi + t* < b\ + Vi, i < n - 1, so 
the original system has at most the same number of customers as the 
modified system before the arrival associated with t' n . Thus, t n > t' n 
cannot occur. Hence t n < t n as claimed. 
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To continue the induction proof for b n , note that (10) and (12) imply 
that di < di for i < n - s. Then, from (11) and (13), we have 

b n = max{£„, d n - s } 5 max{c„, d n - 8 } ^ b' n , 

which completes the proof. □ 

Remark: Our proof of Theorem 8 is closely related not only to Wolff's 
proof, 9 but also to Sonderman's comparison proofs. 20,21 Sonderman was 
concerned with the effect of different service-time distributions instead 
of different queue disciplines. 

We close this section with another result about pure-loss systems. 
With waiting rooms or with general service times it is easy to show 
that the stochastic processes representing the number of customers in 
the system need not be stochastically ordered, but we do get stochastic 
order with exponential distributions and no waiting rooms. 
Theorem 9: In the setting of Theorem 7, if there are no waiting rooms, 
if the service-time distribution is exponential and ifQ(0) = a t Qi(0) + 
Q 2 (0), then 

{Q(t), t > 0} =>., {QAt) + Q 2 (t), t > 0} 

and 

{N{t), t > 0} < s , (iVi(f) + N 2 {t), t > 0}. 

Proof: Here the argument follows Sonderman 7,20,21 and Whitt. 8 As the 
first step in constructing the two systems on the same probability 
space, we let the two systems being compared have identical arrival 
processes; i.e., we let the arrival process to the combined group of Si 
+ s 2 servers be the sum of the two arrival processes to the separate 
groups of Si servers. This not only means that the arrival processes 
have the same joint distributions, but that they have the same sample 
paths. Similarly, we let both systems start off with the same number 
of customers in the system; i.e., given the pair [Qi(0), Q 2 (0)], we let 
Q(0) — Qi(0) + <?2(0). We now show how to construct the departures 
so that 

N(t) < Ni(t) + N 2 (t) (14) 

and 

Q(t) > Qi(t) + Qsit) (15) 

for all t > 0. We generate departures from both systems using a single 
Poisson process with rate (»i + s 2 )fi. Each point in this Poisson process 
corresponds to a potential departure. Suppose the point occurs at time 
L With probability Qi{t)/{si + s 2 ), the point corresponds to a departure 
from both the single group of si servers and the combined group of Si 
+ s 2 servers; with probability Q 2 {t)/(s^ + s 2 ), the point corresponds to 
a departure from both the single group of s 2 servers and the combined 
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group of Si + s 2 servers; with probability [Q(t) — Qi(t) — Q 2 (£)]/(si + 
s 2 ) the point corresponds to a departure from only the combined group 
of Si + s 2 servers; and finally, with probability [si + s 2 — Q(t)]/{si + 
s 2 ), the point corresponds to no departure at all. This can be shown to 
yield the proper distributions for each system; see Sonderman 7 for 
more detail. This also guarantees that there is a departure in the 
combined group of si + s 2 servers whenever there is a departure from 
one of the groups with S] and S2 servers. There also cannot be a 
departure from the combined group alone when Q{t— ) = Qi{t— ) + 
Q2(t—),so inequality (15) is maintained. This means that all departures 
and losses from the combined group of Si + s 2 servers that are not 
matched by corresponding departures or losses from one of the groups 
of si and s 2 servers can be matched with earlier losses from one of the 
groups of s\ and s 2 servers. Mathematical induction on the arrival 
index establishes ( 14) and (15) and formally completes the proof. D 
Remark: For the special case of M/M/s systems, the stochastic order 
in Theorems 6 and 9 can also be established under the conditions in 
the corollary to Theorem 5 using existing comparison theorems for 
continuous- time Markov chains; see Sonderman. 7 However, we know 
of no direct connections between the monotone likelihood-ratio order- 
ings and the sample-path orderings. 

IV. DIFFERENT SERVICE RATES 

In this section we let the service rates in the two separate systems 
be different. One way to extend (2) and (4) occurs when the service 
times are associated with the arrivals. If two independent Poisson 
streams with rates Ai and A 2 and associated service-time c.d.f.'s Fi(x) 
and F 2 (x) are combined, then the resultant stream is a Poisson stream 
with rate Ai + A 2 and associated service-time c.d.f.: 

F{x) = [AiFi(x) + A 2 F 2 U)]/(Ai + A 2 ). 

Of course, when F,(jc) is exponential with mean /i, -1 for each i, F(x) is 
not exponential unless \i\ = \l 2 . However, the blocking probability for 
an M/G/s loss system depends only on the mean service time. Thus 
the loss rate for the combined system is L(si + s 2 , A: + A 2l (Ai + A 2 )/ 
(ai + a 2 )), and a natural extension of (2) to conjecture is 

L(si + s 2 , Ai + A 2 , (Ai + A 2 )/(<2t + a 2 )) < L(si, Ai, /ii) + L(s 2 , A 2 , ji 2 ). 

Unfortunately, this conjectured inequality is not valid. To see this, let 
Xj = \ t ^! = e -1 , X 2 = e, and pi 2 = e 2 ; then Oi = e and a 2 = e _1 . Obviously, 

L(si + s 2 , Aj + A 2 , (Ai + A 2 )/(ai + a 2 )) = (Ai + \ 2 )B(si + s 2 , a x + a 2 ) 

= (1 + e)B(si + s 2 ,e+ e" 1 ) 

^►1 as e -* 0, 
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whereas 

L(si, Ai, Mi) + L(s 2 , A 2 , M2) = Ai5{si, ai) + A 2 B(s 2 , a 2 ) 

- B(si,e) + eB(s 2 ,e- 1 ) 
-h> as e -» 0. 
Consequently, in this case 
L(si + s 2 , Ai + A 2 , (Ai + A 2 )/(ai + a 2 )) > L(«i, A x , jki) + L(s 2 , A 2) ju 2 ) 

for sufficiently small e. The previous measure, rate of customer loss, is 
not the only reasonable way to evaluate system performance in this 
case. For example, one might be interested in the rate of loss of service 
time. (Note that there is no real difference between these measures 
when the mean service times of the systems are identical.) With this 
new measure, the natural extension of (2) to conjecture is: 

° l + ° 2 L(si + s 2 , Ai + A 2) <Ai + A 2 )/(oi + 02)) 
Ai + A 2 

< — L(si, Ai, fii) + — L(s 2 , A 2 , fi 2 ). 

This inequality is in fact always true, since substitution of L(s, A, a) 
= \B(s, X/fi) quickly reduces it to the second version of the inequality 
of Theorem 1. Thus the server occupancy is always increased for the 
combined system. 

Turning to delay systems, we again find examples where sharing can 
be counterproductive. To see this, consider two M/M/l delay sys- 
tems with Ai = 1, /ii = 2, A 2 = e, and a 2 = 2e. Then EQi(co) = EQ 2 {*>) 
= p/(l — p) = 1, but EQ{oo) can be shown to be of order e -1 as e -* 0: 
Consider the interval following a low-intensity arrival. With probability 
A 2 /(A 2 + /x 2 ) = V6, a second low-intensity arrival occurs before the first 
departs. Then there follows an exponentially distributed interval of 
mean length l/4e during which the combined system fills up with high- 
intensity customers. In computing the average number of customers 
in the system, we get a term of order e -2 (the total area in the plot of 
the number of customers in the system versus time, starting from the 
moment the second low-intensity customer arrives and ending when 
one of the two low-intensity customers departs), divided by a term of 
order e"\ In other words, with the mean steady-state delays held fixed 
in the two separate systems, the mean steady- state delay in the 
combined system can be arbitrarily large. 

Note that the combined system can be modeled as an M/G/si+s 2 
delay system where the service-time distribution is the mixture of two 
exponential distributions, but in contrast to the case of loss systems 
the mean delay does not depend only on the mean of the service-time 
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distribution. Hence, the appropriate generalization of (4) involves a 
system which is not M/M/s. 

Another possible extension for fi\ ¥> fi2 occurs when the service-time 
distributions are associated with the servers. Here the combined 
system is not M/M/s because there are heterogeneous servers, so there 
are no equations similar to (2) and (4). In this case, it can be shown 
that with exponential service-time distributions and no waiting rooms, 
resource sharing is always better if customers always are sent to the 
fastest available server. In particular, as in Theorems 6 to 8, it can be 
shown for any single system that assigning customers to the fastest 
available server produces fewer losses than any other rule, where by 
"fewer losses" we mean in the sample-path ordering of Section III. One 
other rule, corresponding to the two separate systems, is to assign the 
customer only to servers associated with their original separate arrival 
streams. 

When we focus on delay systems with heterogeneous servers, it is 
easy to give counterexamples showing that resource sharing can again 
be counterproductive. Related literature on the assignment of cus- 
tomers to heterogeneous servers appears in Winston, 22 " 24 Smith, 25 and 
references therein. 

This section shows that, with unequal service-time distributions, 
resource sharing can be counterproductive. However, with unequal 
service-time distributions, much depends on the criterion of system 
performance. Also, it should be noted that such counterexamples have 
been observed before; others have discovered that infrequent "bad" 
customers can affect adversely a large number of "good" customers. 
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APPENDIX 

Here we give two results that are due to others. First, we present Paul 
Burke's proof that B(ts, ta) is strictly decreasing in t for t > 0. This 
result implies Theorem 1 when Ai/si = a 2 /«2, because then 
B(si + s 2f a\ + a 2 ) = B(tSi, tat) for some t > 1, so B(s,, a*) > 
B(si + s 2 , a\ + a 2 ) for each i and 

AiB(si,ai) A 2 B(s 2 ,a 2 ) 
i ^ i — + — r- — r — ^ B(si + s 2 , cti + a 2 ), 

A] + A2 Ai t A 2 

which is equivalent to (2). 

To see that B(ts t ta) is strictly decreasing in t t first recall the 
following equation relating two different expressions for the tail of the 
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gamma distribution: 

k=n f" „n -x 

e -« T a */A!- — i-rf*- 



Then note that 



1 



e-V d* 



B<fc>, to) e- ta {ta) ts 



-I. 



e-^-^'tx/to)'" d* 

e _jr U + [x/ta])'" dx; 
Jo 

also see Theorem 3 of Jagerman. 3 Finally, (1 + [x/ta])' 8 is strictly 
increasing in t. 

Second, Herbert Shulman has shown that Paul Burke's result and 
the convexity of B{s, a) in s for s > 1 imply a version of Theorem 1. 
Such convexity has frequently been conjectured but has been proved 
only for lattices of points with unit spacing, see Messerli 26 and refer- 
ences therein. These versions of convexity are not strong enough to 
make the following valid even when si and s 2 are integers; however, 
general convexity would establish the proof for all real numbers si and 
$ 2 > 1, a more general mathematical result. We reproduce Shulman's 
argument here: 

cti ~/°i + a * \ 

B(s y + s 2 , a r + a 2 ) < B[ Si , a, + a 2 

a, + ai \ a\ j 

a-i n fa } + a a . \ 

B[ s 2 , ai + az\ 

ai + a 2 \ ai J 



+ 



B(gi,ai)+ ° 2 B(s 2t a 2 ). 



a\ + a-i a\ + a-i 
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